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ABSTRACT 

Power quality in AC power system depends on the harmonic 
content of the waveform. To accurately predict the content of 
harmonic components in a particular system configuration, an 
engineer needs the analytical tools to predict the behavior of the 
various harmonic components. This paper describes a harmonic 
power flow program applicable to space power systems with 
sources of harmonic distortion. The algorithm is a modification of 
EPRI’s HARMFLO program which assumes a three phase, 
balanced, AC system with loads of harmonic distortion. The 
modified power flow program can be used with single phase, AC 
systems. 

Early results indicate that the required modifications and the 
models developed are quite adequate for the analysis of a 20 kHz 
testbed built by General Dynamics Corporation (GDC). This is 
demonstrated by the acceptable correlation of the present results 
with published data. Although the results are not exact, the 
discrepancies are relatively small. 


INTRODUCTION 

As the size of spacecraft has grown over the years, so has the 
demand for electrical energy consumed by these vehicles. In the 
first manned spacecraft, Mercury, the power supply was a low 
capacity, 24 VDC system. As the plans for the vehicle needed to 
reach Mars indicate, there must be more capability to support 
more people and more equipment for a much longer period of 
time than any of the previous projects have in the past. This 
capability is greatly enhanced by an AC system. The proposed 
high frequency 20 kHz system [1] has the advantages of high 
efficiency and reliability, minimum overall system weight, and easy 
system growth. Although an AC system can solve the power 
demand problem, it can cause problems which must be 
addressed to assure an efficient, reliable power supply. 

A potential problem with an AC system is harmonic resonance. 
Since the proposed source of the 20 kHz, AC wave is an inverter 
and the system will supply various nonlinear loads among others, 
harmonic resonance is a possible occurrence. Computer software 
is required to study this phenomenon so that proper filtering can 
be designed for its elimination. 


*NASA/ASEE Summer faculty fellow at NASA Lewis Research Center 


Voltage and current resonance in a power system can cause 
many serious problems. These include: 

• insulation failure due to overvoltage, 

• false instrumentation data due to multiple zero crossings of 
the waveform, 

• current overload in lines and other components, 

• rf noise, and 

• malfunction of electrical equipment due to high frequency 
signals in the power supplied. 

Other problems created by low power quality in the power supply 
system may exist. However, the above problems are intended to 
demonstrate the reason for concern over harmonic resonance. 

A reasonable method for identifying and evaluating resonance 
problems in a proposed power system is to use a harmonic power 
flow program such as HARMFLO developed by the Electric Power 
Research Institute (EPRI). HARMFLO can predict steady-state 
resonance conditions in three phase, balanced power systems at 
60 Hz. 

This paper discusses a modified harmonic power flow algorithm 
for use with mainly spacecraft power systems. Unlike its terrestrial 
counterpart used for a balanced, three phase, 60 Hz system with 
loads of harmonic distortion, the algorithm described here is 
applicable to a high frequency, single phase, AC system with a 
source of harmonic distortion. 


THE HARMONIC POWER FLOW ALGORITHM 

The HARMFLO algorithm, developed by Dr. G. T. Heydt [2,3] of 
Purdue University, uses a modified version of the standard 
Newton-Raphson formulation for a power flow solution. The 
modification involves the alteration of the voltage vector to include 
all of the harmonic components to be examined not just the non- 
triple odd harmonics typically found in three phase balanced 
power systems. Also, non-existent models of components such as 
the Mapham inverter and a single phase rectifier have been 
added. The voltage vector becomes 

[ y<D, ... v<». ... V<">. $> ] r 

where V * 1 * is the set of i th harmonic bus voltages and 
represents the firing and commutation angle of any three phase 
rectifiers. The formulation of the harmonic power flow also 
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requires the modification of the scheduled quantities to include 
those produced by the inverters. 

A /(I ) = j(l )*ch*dubd _ i(l)c*cui«*d 


where I • ch#du,#d js the scheduled i th harmonic current at each 
bus, and I ^ c * ,cu, * ud is the calculated i th harmonic current at 
each bus. 

The process also requires the modification of the Jacobian matrix 
to 
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Thus, this modified HARMFLO differs from others which use the 
fundamental power and harmonic current responses of the 
nonlinear devices to solve the voltage levels within the power 
system [4]. Also, it differs from the use of ATP (Alternate 
Transients Program) version of the Electromagnetic Transient 
Program (EMTP) to perform harmonic analysis in which nonlinear 
loads are represented by harmonic current injections at desired 
nodes within the power system [5]. 


MODIFICATION TO THE HARMONIC POWER FLOW ALGORITHM 

The original harmonic power flow algorithm, HARMFLO, cannot 
be used directly for the study of spacecraft power systems for two 
major reasons. The first reason is that the algorithm, as 
developed, is designed strictly for three phase, balanced power 
systems rather than single phase systems. Secondly, only loads, 
but not the power sources, can be sources of harmonic distortion. 
Both of these assumptions are adequate for terrestrial power 
systems. However, spacecraft power systems are designed to be 
single phase with solid state DC-AC inverters. These potential 
harmonic sources imply the need for modifications to the basic 
algorithm of HARMFLO. 

The modifications needed to handle all harmonic components 
are relatively simple. They comprise redimensioning of the 
matrices, and adjustment of pointers in the algorithm to assure the 
exclusion of undesired harmonics from the analysis. 

The harmonic power flow algorithm assumes that all AC power 
is supplied by balanced, three phase sinusoidal sources. 
However, a potential power source for spacecraft can be an AC 
source which embodies DC-AC solid state inverters. The most 
commonly cited AC source to employ the resonant type of 
operation is the Mapham inverter [6] shown in Figure 1. The 
Mapham inverter, popular due to its simplistic, light weight design 



Figure 1 - Circuit diagram of a Mapham Inverter. 


and good regulation, employs an L-C series combination to 
produce a dampened AC waveform. The SCR’s are used to 
reverse the current on every half cycle to keep the inverter 
producing AC power. The resulting output of this type of inverter 
is shown in Figure 2. The discontinuities in the waveform indicate 
that this power source can cause harmonic distortion. This means 
that, to use HARMFLO in this application, the Mapham inverter 
must be modelled in the harmonic power flow algorithm. The 
required modifications for modelling such a power source are 
relatively more complex. 
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REGIONS FOR MODELLING 



Figure 4 - Capacitor Current, I ^ Inductor Currents, I L1 and 
I ^ with Indicated Regions for Modelling. 


MODELLING THE MAPHAM INVERTER 

The modelling of this type of device requires the computation 
of the scheduled magnitude and phase angle of each harmonic 
current component at each iteration. The first step in this process 
is to find the resonant frequency of the RLC circuit formed by the 
Mapham inverter and the power system. The power system 
equivalent determined at the terminals of the Mapham inverter is 
added to the inverter model as shown in Figure 3. Ordinary circuit 



Frequency of System. 


analysis techniques, such as Laplace transforms, can be used to 
find r r and u r , which are the time constant and frequency, 
respectively, of the natural response of the circuit. 

The values of r r and u r can then be used to find an analytical 
representation of one period of the current waveform. This is 
accomplished by dividing a period of the waveform into regions 
defined by the physical characteristics of the waveform. Since the 
Mapham inverter is operated by firing the SCR's for the next half 
cycle before the conducting SCR's have commutated, the current 
waveform is defined by five distinct regions as shown in Figure 4. 
The resulting analytical expressions for each of the regions shown 
in Figure 4 are: 


SCR 2 and SCR 3 are conducting and SCR , and SCR 4 
are off, 

lc,(t) = -l 0 e \ 2 'sin 


for 0 it< hZJl ; 
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SCR v SCR 2 , SCR 3 , and SCR 4 are conducting, 
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5. SCR 2 and SCR 3 are conducting and SCR 1 and SCR 4 
are off, 
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where a is the damping coefficient; r $ is the period of the 
system frequency. The value of 1 0 is determined at each iteration. 
However, since the system configuration is fixed, r t , r r , and « t 
do not change in value during solution. 

Fourier transforms are then applied to the five regions to 
determine the magnitude coefficient and phase angle for each 
harmonic component. These values are inserted into the 
scheduled value matrix for the appropriate busses. The algorithm 
is then allowed to proceed as it normally would. 


TEST RESULTS 

To demonstrate the theories presented in this paper, the 
General Dynamics Corporation 20 kHz testbed was used as a 
source of comparison [7]. The testbed was configured to have a 
Mapham inverter supply a constant power AC load through a 50 
meter transmission line. The Mapham inverter has the following 
characteristics: 


Operating Voltage: 

440 

V. 

System Frequency: 

20 

kHz, 

L,= 

16 

MH, 

C,= 

1.71 

mF, 

C,: 

2.0 

H F. 


The transmission line has the following characteristics: 

R: 1 .043 m Q/meter, 

L 0.027 ilH/ meter, 

C: 0.003 jiF/meter. 

The characteristics of the load are: 

Power: 2-10 kVA, 

p.f.: 80%, lagging. 


The resulting Total Harmonic Distortion (THD) of the bus voltage, 
defined as 


THD 


E 
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v y) 1 


was obtained from the harmonic power flow studies and these 
results are compared with actual test data from the testbed. The 
results are summarized in Rgure 5. A resonable correlation is 
evident between the experimental results and those obtained from 
the modified HARMFLO. 



Figure 5 - Comparison of Results of Bus Voltage THD 
versus KVA Loading from Harmonic Power Row 
Program and Test Data from the GDC 20 kHz 
Testbed. Current THD versus KVA Loading is 
also included. 


CONCLUSIONS 

The simulation of harmonic behavior in large spacecraft power 
systems can now be studied and problems associated with this 
behavior can be corrected since a major analytical tool needed for 
this work is now available. 

Rne tuning of the system models and their values, via 
experimental verification, are the items which must be further 
evaluated so that accurate conclusions can be drawn on the 
complete validity of the source code. The completion of this 
evaluation work is scheduled for 1992. 

The models used are flexible enough to allow for future 
changes in system design and operation. The modelling has also 
been incorporated into HARMFLO with little degradation in the 
convergence characteristics of the Newton-Raphson method. 

The modified HARMFLO can be used for the determination 
and/or prediction of system interactions such as harmonic 
resonance and electromagnetic interference (EMI). A future paper 
will address these issues. 
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